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We discuss a robust projection method for the extraction of excited-state masses of the nucleon 
from a matrix of correlation functions. To illustrate the algorithm in practice, we present results for 
the positive parity excited states of the nucleon in quenched QCD. Using eigenvectors obtained via 
the variational method, we construct an eigenstate-projected correlation function amenable to stan- 
dard analysis techniques. The method displays its utility when comparing results from the fit of the 
projected correlation function with those obtained from the eigenvalues of the variational method. 
Standard nucleon interpolators are considered, with 2x2 and 3x3 correlation matrix analyses pre- 
sented using various combinations of source-smeared, sink-smeared and smeared-smeared correlation 
functions. Using these new robust methods, we observe a systematic dependency of the extracted 
nucleon excited-state masses on source- and sink-smearing levels. To the best of our knowledge, this 
is the first clear indication that a correlation matrix of standard nucleon interpolators is insufficient 
to isolate the eigenstates of QCD. 
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I. INTRODUCTION 

Lattice QCD provides a non-perturbative tool to ex- 
plore many properties of hadrons from first principles. 
In the case of the hadron mass spectrum there are well 
developed methods to compute the mass spectra. How- 
ever, while the extraction of the ground states of the 
hadron spectrum is a well understood problem, and has 
provided impressive agreement with experimental results 
[H, the excited states still prove a significant challenge. 
The Euclidean-time correlation function provides access 
to a tower of states since it is a sum of decaying expo- 
nentials with the masses of the states in the exponents. 
The ground state mass, being the lowest energy state 
and thereby having the slowest decay rate, is obtained by 
analysis of the large time behaviour of this function. The 
excited states, however, belong to the sub-leading expo- 
nentials of the two-point correlation function. Extracting 
excited states masses from these sub-leading exponents 
is difficult as the correlation functions decay quickly and 
the signal to noise ratio deteriorates more rapidly. 

One of the long-standing puzzles in hadron spec- 
troscopy has been the low mass of the first positive par- 
ity, = i^, excitation of the nucleon, known as the 
Roper resonance A^*(1440 MeV). In constituent or va- 
lence quark models with harmonic oscillator potentials, 
the lowest-lying odd parity states naturally occurs be- 



low the N 



1 + 



state (with principal quantum number 

er resonance 



iV = 2) 0, [f whereas, in nature the Rop^ 

is almost 100 MeV below the N = ^"(1535 MeV) state. 
Similar difficulties in the level orderings appear for the 
= 1"^ A* (1600) and i"^S*(1690) resonances, which 
have led to the speculation that the Roper resonance may 
be more appropriately viewed as a hybrid baryon state 
with explicitly excited glue field configurations [1, Q 



or as a breathing mode of the ground state [6|] or states 
which can be described in terms of meson-baryon dynam- 
ics alone Q. The first detailed analysis of the positive 
parity excitation of nucleon was performed in Ref. Q 
using Wilson fermions and an operator product expan- 
sion spectral ansatz. Since then several attempts have 
been made to address these issues in the lattice frame- 
work [i, [m, [H, E [B, m, Hi, but in many cases 
no potential identification of the Roper state has been 
made lol ficj. Il,|3[l3l- Recently however, in the anal- 
ysis of [ij, [H, [13], a low- lying Roper state has been 
observed by using advanced fitting techniques [H, [ll| 
based on Bayesian priors. Significant finite volume effects 

on the first pos itive parity state have been observed 
in Refs. [la.l20l.l2ii using the Maximum Entropy Method 
p^. [19I. [22l. |23|. [241 . Here, we use another state-of-the-art 
approach, namely 'the variational method' [25l. [26l. [27l| . 
which is based on the correlation matrix analysis and has 



trix analysis and has 



been used quite extensively in Refs . . . . 

M, M, l32, 33, M, 35, M, M M, m,m m m, El with 

the first analysis of the nucleon performed by Sasaki et al 
. Though the ground state mass of the nucleon has 
been described successfully, an unambiguous determina- 
tion of the Roper state has not been successful to date 
with this method, though significant amounts of research 
have been carried out in Ref. [Ill, HI] , the CSSM Lattice 
Collaboration [H, [Mi , the BGR [U, [^, [H, [H, [H 
collaboration and in Refs. [ZH. |4^. 

In this paper, we discuss an analysis method to extract 
masses of the nucleon from the correlation functions using 
a variational analysis. Employing standard interpolating 
operators xi, X2 and X4, we discuss the method for 2 x 2 
and 3x3 correlation matrices with the point and a range 
of sweeps of Gaussian smearing [3] at the source, sink 
and at both source and sink. This analysis shows for the 
first time (and despite the fact that people have been 
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using source smeared correlation functions for quite long 

time [il, M, M, M, M, M, M, M, M, M, S, [la, S, m ) 

that, unexpectedly, the excited states of the nucleon are 
smearing dependent. This analysis indicates that signif- 
icant caution should be taken when employing a partic- 
ular level of smearing. To ensure the maximal indepen- 
dence of our results on human input and minimal errors, 
we construct a 'robot' algorithm, governed by defined fit- 
ting criteria, that automatically performs a standardised 
fitting procedure. We present here results from this al- 
gorithm and also those obtained from the eigenvalues to 
provide confidence in the extraction of the nucleon mass 
spectrum. 

This paper is arranged as follows: Section ITTl contains 
the general description of the extraction of masses with 
the introduction of different nucleon interpolating fields. 
The lattice details are in Section Hill the analysis method 
is presented in Section HVl and conclusions are presented 
in Section IVl 

II. MASS OF HADRONS 

The masses of hadrons are extracted from two-point 
correlation functions using operators chosen to have over- 
lap with desired states. Let us consider a baryon state 
B of spin half, if we suppress Dirac indices a two point 
function can be written as, 

G,,{t,p) = X^e-'^^-~(r!|T{x.(x)x,(0)}|17). (1) 

X 



The operator Xj(0) creates states from the vacuum at 
space-time point and, following the evolution of the 
states in time t, the states are destroyed by the operator 
Xi(a;) at point x. T stands for the time ordered prod- 
uct of operators. Having a complete set of momentum 
eigenstates requires that, 

^ \B,p\s){B,p',s\^I, (2) 

where B can include multi-particle states. The substitu- 
tion of Eq.jH) into the Eq.jT]) yields, 

G,,{t,p)^J2 E e-'P■'{n\x^{x)\B,p',s){B,P,s\x,mn). 

S B,p',s 

(3) 

We can express the operator as 

X^{x) = e'^^■-x.(0)e-'^■^ (4) 

where, P'^ = P = {H, P) and P is the momentum op- 
erator whose eigenvalue is the total momentum of the 
system. Eq.([ni) can now be written as. 



J 



G.j{t,p)^Y. E e-^'^■'{^\e'''^X^me-'''''\B,p',s){B,P,s\xJ{Qm 

= E E e-^''''e-^'■'^^-^^'H^\x^{m.P' :S){B,p' ,s\xAm)- (5) 



As we move from Minkowski space to Euclidean space, states by 
the time t — > —it and the above equation then can be 
written as, 



(17|x(0)|i3+,p,s) = Ab+ (7) 



E 



G^,it,p)= X] e-^-\p-,(l}|x.(0)|i?,p',s)(S,p*',sfc(0)|f7) ' 

B,p',s 

. X;Ee-'(n|x,(0)|B,p..)(i,,S,|x,(0)|«). = A,.,/^««.(p». (8) 

(6) 

For the negative parity states the definition is 

The overlap of the interpolating fields x(0) and x(0) with 
positive and negative parity baryon states \B^) can be 

parametrized by a complex quantity called the coupling {^\x{0)\B~ ,p, s) = A^- J -^—Ib'^B- (Pi s), (9) 

strength, Ab±, which can be defined for positive parity V 
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(i?-,p,s|x(0)|17) = -AB-. 



-UB-ip,s)"f5. (10) 



Here, Xb± and \b± are the couplings of the interpolat- 
ing functions at the sink and the source respectively and 
Mb± is the mass of the state . Eb± is the energy of 

the state B^, where Eb± = \J A/^± + p^, and ub±{p, s) 

and ub± {p, s) arc the Dirac spinors, 



(11) 



Thus, Eq.Q contains a projection operator r± = 
'^sU^±{p, s)u'^±{p, s), through which the contributions 
to the even and odd parity states from the correlation 
function can be obtained. For positive parity, this can 
be expressed as, 



8 

and for the negative parity. 



7.p + Mb+ 



2E 



B+ 



(12) 



75 



^ <_ (p, s)u%^ ip, .) 75 = (13) 



2Ef 



By substituting the above Eqs. for the positive and neg- 
ative parity states in Eq.Q we obtain. 



,tl-PB^ 



M 



B + 



B + 

E 

B 



2Eb+ 

E„_t ~l-PB- +Mb- 

2Eb- 



(14) 



At momentum p = 0, Eb± = Mb±, and a parity projec- 
tion operator r± can be introduced, 



r± = ^(i±7o) 



(15) 



We can isolate the masses of the even and odd parity 
states by taking the trace of G with the operators r+ and 
r_. The positive parity state propagates through the 
(1,1) and (2,2) elements of the Dirac matrix, whereas, 
negative parity state propagates through the (3,3) and 
(4,4) elements. 

The correlation function for positive and negative par- 
ity states can then be written as. 



G±(t,0)=Tr,p[r±ay-(t,0)] 



B± 



(16) 



The correlation function contains a superposition of 
states. The mass of the lowest state, Mq± can be ex- 
tracted at large t where the contributions from all other 
states are suppressed. 



G±(t,0) 



(17) 



The source smearing |44| technique is applied to in- 
crease the overlap of the interpolators with the lower ly- 
ing states. A fixed boundary condition in the time di- 
rection is applied for the fermions by setting Ut{x^ Nt) = 
OVx in the hopping terms of the fermion action with 
periodic boundary conditions imposed in the spatial di- 
rections. Gauge invariant Gaussian smearing [44j in the 
spatial dimensions is applied through an iterative pro- 
cess. The smearing procedure is: 



i{}i{x,t) = YF{x,x')ipi-i{x' ,t), 



(18) 



where. 



F{x,x') = (1 - a)Sx,x' + - ^[?7^(a;)4',2;-H/i 

+ Ulix~fi)S,,^,^f,], (19) 

where the parameter a = 0.7 is used in our calcu- 
lation. After repeating the procedures iVsm times on a 
point source the resulting smeared fermion field is. 



i;N,Jx,t) = ^F^-(x,x')1^o(a;^^)■ 



(20) 



The extraction of the ground state mass is done 
straightforwardly. However access to the excited state 
masses requires additional effort. Here we consider the 
variational method [1^, [2^, [l^ . The variational method 
requires the cross correlation of operators so that the op- 
erator space can be diagonalised and the excited state 
masses extracted from the exponential nature of the di- 
agonalised basis. To access N states of the spectrum, 
one requires a minimum of N interpolators. Tradition- 
ally, only a few interpolators are considered providing 
access to a small number of states of the desired channel. 

The parity projected two point correlation function 
matrix for p ~ can be written as, 

G,,(t) = (X]Tr,p{r±(r!|x.(x)x,(0)|f^)}) (21) 

X 

N-1 

= Y K^e-^'-'- (22) 

Here, Af and A" are the couplings of interpolators Xi 
and Xj at the sink and source respectively to eigenstates 
a = 0, • • • , (iV — 1). rua is the mass of the state a. The N 
interpolators have the same quantum numbers and pro- 
vide an A^-dimensional basis upon which to describe the 
states. Using this basis we aim to construct N indepen- 
dent interpolating source and sink fields which isolate A'^ 
baryon states \Ba), i.e. 



N 

E 



Ui Xi 



(23) 
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N 



r=Y.v^X^. (24) 

i=l 



such that, 



{Bp,p,s\4>'^\n) = 5^fiz''u{a,p,s), (25) 
{n\<j>''\Bp,p,s) ^ 5^pz''u(a,p,s), (26) 

where and z" are the coupUng strengths of 0" and 0" 
to the state 1-8^). Consider a real eigenvector which 
operates on the correlation matrix Gij(t) from right, one 
can obtain H^ . 



G,:,(t)^," = (5^Tr,p{r±(r!|x.X,|r!)})i 



Q -OL ^—jrictt 



(27) 



For notational convenience, in the remainder of the dis- 
cussion the repeated indices i, j, fc arc to be understood as 
being summed over, whereas, a, which stands for a par- 
ticular state, is not. Since the only t dependence comes 
from the exponential term, we can write a recurrence re- 
lation at time {i. + At) as, 



Gy (t + At)w" ^ e-™°^'Gy(tX 



(28) 



If wc multiply the above equation by [Gy (t)] ^ from the 
left wc get, 

[{G[t))-^G{t + At)]u" = e-^-^'u" 

= c"u". (29) 

This is an eigenvalue equation for eigenvector with 
eigenvalue c" = e"™""^*. We can also solve the left eigen- 
value equation to recover the v°' eigenvector, 



<Gy(i + At) = e-"°^*<Gy(t). (30) 



Similarly, 



z;"[G(t + At)(G(t))-i] 



(31) 



The vectors and diagonalize the correlation matrix 
at time t and t + At making the projected correlation 
matrix. 



(32) 



The parity projected, eigenstate projected correlator, 
v°'Gfj{t)u'j = G% is then used to obtain masses of differ- 
ent states. We construct the effective mass 



Me"ff(t) = In 



G^(t,0) 
G5(t + 1,0), 



Ml. 



(33) 



and apply standard analysis techniques as described in 
the following. 



III. LATTICE DETAILS 

We use an ensemble of 200 quenched configurations 
with a lattice volume of 16'^ x 32. Gauge field configu- 
rations are generated by using the DBW2 gauge action 
[13, [EH and an C'(a)-improved FLIC fermion action [s^l 
is used to generate quark propagators. This action has 
excellent scaling properties and provides near continuum 
results at finite lattice spacing [53| . The lattice spacing is 
a = 0.1273 fm, as determined by the static quark poten- 
tial, with the scale set with the Sommer scale, Tq — 0.49 
fm |54l | . In the irrelevant operators of the fermion action 
we apply four sweeps of stout-link smearing to the gauge 
links to reduce the coupling with the high frequency 
modes of the theory . We use the same method as in 
Ref. [4^ to determine fixed boundary effects, and the ef- 
fects are significant only after time slice 25 in the present 
analysis. Beside point operators, various sweeps (1, 3, 7, 
12, 16, 26, 35, 48, 65 sweeps corresponding to rms radii, 
in lattice units, 0.6897, 1.0459, 1.5831, 2.0639, 2.3792, 
3.0284, 3.5237, 4.1868, 5.0067) of gauge invariant Gaus- 
sian smearing [IJ are applied at the source (at t = 4) 
and at the sink. This is to ensure a variety of overlaps of 
the interpolators with the lower-lying states. The anal- 
ysis is performed on four different quark masses provid- 
ing pion masses {0.797,0.729,0.641,0.541} GeV. 
The error analysis is performed using a second-order sin- 
gle elimination jackknife method, where the x^/dof is 
obtained via a covariance matrix analysis method. We 
discuss our fitting method in the next section. 

The nucleon interpolators we consider in this analysis 
are, 

Xi(x) = e"^^(«^"(x)G75d'(a:))w^(x) (34) 
X2{x) = e^''^{u^''{x)Cd\x))^5u'{x) (35) 
X4(x) = e''''%u^''{x)Cj5l4dHx))u%x) (36) 

The xi and X2 interpolators arc used in Refs. [sl. [ill. [56j . 
The X4 interpolator is considered as the time component 
of the X3 interpolator used in Refs.jH, [U, [l^]. We use 
the Dirac representation of the gamma matrices in our 
analysis. 



IV. ANALYSIS AND DISCUSSION 

A. Correlation Matrix Analysis 

Wc begin this section by outlining the fitting method 
used in this paper. The method is based on the maximi- 
sation of the Euclidean-time fit window and the minimi- 
sation of the value of the x^/dof commencing at the ear- 
liest time possible (s^ • At larger times the lighter states 
dominate the correlation function. However, the signal 
to noise ratio decreases rapidly, forcing increased errors 
in the results. Additionally, ignoring smaller time slices 
may lead to eliminating important information about ex- 
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FIG. 1: (Color online). Mass of the nucleon (A''^ ) from the projected correlation function as shown in Eg 1321 fleft) and from 
the eigenvalue (right) for a 2 x 2 correlation matrix of xi, X2 interpolators. The figure corresponds to a pion mass of 797 
MeV (heaviest) and for the point source to point sink correlation functions. Each pair of ground and excited states masses 
correspond to the diagonalization of the correlation matrix for each set of variational parameters tstart (shown in major tick 
marks) and At (shown in minor tick marks). Here, the time t as shown in Eas. (|29l and 13 If) is called tstart. 
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FIG. 2: (Color online). As in Fig.[T] but for source smeared to point sink correlation functions with the number of Gaussian 
smearing sweeps set to 7, which corresponds to a smearing radius of 1.5831 in lattice units. 



cited states included in the two-point correlation func- 
tion [sll. Data at larger time slices for excited states 
extracted using the variational method can be contam- 
inated by residual contributions from the lower lying 
states, resulting in lower values of excited states mass. 
In consideration of these points, in this analysis we use 
a preference towards earlier times which have a high sig- 
nal to noise ratio and are therefore heavily constraining 
in the fit procedure. This allows us to isolate the energy 
level from higher state contamination through the x^/dof 
and simultaneously control the errors potentially intro- 
duced at higher times (where the signal to noise ratio is 
lower) by contamination from lower-lying states in the 
varitiational analysis procedure. Hence, the inclusion of 



these early times minimizes the error in the results while 
an acceptable value of x^/dof is maintained. 
We perform the analysis for 2x2 and 3x3 correlation ma- 
trices via the variational method. For the quark masses 
considered herein, the interpolator xi has better overlap 
with the lower-energy states 113 a nd strongly couples 
to the nucleon ground state [111 . |32| | , whereas the inter- 
polator X2 does not have good overlap with the nucleon 
ground state and couples to the higher energy state(s) 
[1, [nl, in, . We have found that the other interpola- 
tor, X4, is very similar to xi and also couples strongly to 
the ground state (s^ . which suggests that the xi ^^nd Xi 
operators are somewhat linearly dependent to each other. 
We call the start time, t, of the variational analysis tstart- 
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The diagonalization is accomplished for different values 
of tstart with a few values of At, here At =1 — 5 for each 

^start ■ 

We now consider fits to the parity and eigenstate pro- 
jected effective mass of Eg. ([55]) which follows from the 
eigenvector analysis. In this presentation, we consider 
the cut-off value for an acceptable value of x^/dof as 
1.30. Firstly, we try to fit the effective mass from two 
time slices after the source to the largest possible time 
of the correlation function i.e., from time slice 6 (since 
the source is at 4) to time slice 25 (after which the fixed 
boundary effects are significant). We call the lower time 
tmiii and the larger time <max- If an acceptable fit is not 
obtained then we keep tmin fixed and decrease tmax and 
reattempt a fit. If this is also unsuccessful then we it- 
erate the same process until we reach a time near tnmi 
while maintaining a minimum fit window size. At this 
point, if an acceptable fit (as dictated by the x^/dof) 
is still not obtained then we increase imin by one time 
slice and try to fit the new window tmin to tmax- This 
process repeats until a fit is obtained. The minimum fit 
window size we consider for the ground state is 5 time 
slices in the effective mass which corresponds to 6 time 
slices in G{t). For the excited states, the minimum win- 
dow size of 3 is considered corresponding to 4 time slices 
in G{t). This provides a balance in providing evidence of 
an eigenstate while avoiding residual contaminations of 
lower-lying states. 

Figs. [T] and [5] present the ground and excited states of 
the nucleon for a 2 x 2 correlation matrix with xi and 
X2 interpolators for point and smeared source correlation 
functions respectively. The point-like correlation func- 
tion is a difhcult correlator to extract a mass from, as it 
admits strong overlap with excited states. Nevertheless, 
we consider the point correlation function as a challenge 
in this analysis. 

The left and right figures show the mass from the pro- 
jected correlation function and the mass from the eigen- 
value respectively. Each point corresponds to the diago- 
nalization of the matrix for each set of variational param- 
eters istart and At. The mass coming from eigenvalues 
are intrinsic to the variational analysis since they come 
directly from the diagonalization of the matrix, while the 
mass from the projected correlation function comes from 
the 'robot' algorithm described above. It is interesting 
to note that masses from the projected correlation func- 
tions are almost independent of the variational parame- 
ters. Fig. [1] also shows that mass can also be extracted 
reliably from the point-to-point correlation function. The 
behaviour of the eigenvalues at lower igtart and At reflects 
the contamination of higher excited states. Although it 
can be difficult to extract a mass directly from the eigen- 
values, it is relatively easy to expose a mass in the pro- 
jected correlation function. 

Figs. [3] and [4] show the eigenvectors for the diagonali- 
sation of correlation matrices for the point (for Fig. [T|) 
and source-smeared (for Fig. [5]), correlation functions re- 
spectively. Eigenvectors are normalized for each set of 



variational parameters to unit length. It is interesting to 
note that the eigenvectors do not show a strong sensi- 
tivity to excited-state contamination. As with the mass 
from the eigenvalue, at larger tstart, the eigenvectors are 
also dominated by errors. Eigenvectors in Figs. [3] and 
13] also indicate that as the xi and X2 interpolators are 
much orthogonal to each other [l^ , the xi interpolator 
has little influence over the excited state and the X2 inter- 
polator also contributes very little to the nucleon ground 
state. 

Figs. inland [S] show the ground and excited state mass 
of the nucleon (as in Figs. [T] and [2|) for xi and X4 in- 
terpolators. The excited state mass for the point cor- 
relation function. Fig. [5l starts a little from below and 
settles down after a few values of tstart- It should be 
noted that similar effects also persist in the 3x3 corre- 
lation matrix analysis where both xi and X4 operators 
are present. However, this behaviour diminishes with the 
level of smearing as shown in Fig. [51 
The correlation matrix analysis of xi and X4 interpola- 
tors provide nontrivial mixing as illustrated in Figs. [3 
and [51 As discussed earlier, xi and X4 are very similar 
and there is little to separate them. This is illustrated in 
Fig. [71 which shows some drift in the eigenvectors for the 
ground state, but with little variance in the ground-state 
mass. 

Fig. [31 presents eigenvectors of the 3x3 correlation ma- 
trix of xi, X2 and X4 interpolators for the point-to-point 
correlation functions. As xi and X4 interpolators largely 
couple to the ground state, the left figure in Fig. [51 shows 
the higher contributions to the ground state come from 
the Xi and X4 interpolators, while the first excited state 
(middle figure) is completely dominated by X2 interpola- 
tor and in the second excited state (right figure) contri- 
butions from all the interpolators are distributed. 
To select a single mass from a series of tstart and At, a 
value of istart = 8 is preferred, and At > 4 (if possible). 
We prefer to avoid the value of At = 1, as it appears 
that At = 1 is, in a few cases, more prone to fluctuations 
than larger values. If a mass is not obtained for these pa- 
rameters (this is the case for the lighter quark masses), 
then we decrease At by one time slice and try to obtain 
a diagonalisation. This procedure is repeated until the 
value tstart + At = 10 is reached. If the diagonalisation is 
still not obtained then we decrease tstart and repeat the 
same procedure. We note that our experience is in ac- 
cord with that of Ref. [13, [Mi • In practice we emphasise 
the importance of keeping tgtart + At large [6l| . 



B. Smearing Dependency of Excited States 

Now we discuss the smearing dependency we have ob- 
served in the masses of the excited states. In Figs. [H 
and [21 a careful comparison of the masses for point- and 
source-smeared correlation functions reveals that the ex- 
cited state mass for the source-smeared case is lower in 
value than for the point source. In contrast, the ground 
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FIG. 3; (Color online). Eigenvector values for vf (left) and uf (right), as shown in Eq. (|32|) . for the correlation matrix analysis 
of Fig. [T] The superscript a stands for the eigenstates while the subscript i represents the interpolators. Here, for the 2x2 
correlation matrix, a = 1, 2 and i = 1, 2. 
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12 



state masses agree within one standard deviation for al- 
most all sets of variational parameters. 

Here we extend the analysis for various amounts of 
source-smearing sweeps in the correlation functions. Our 
new robust analysis techniques reveal that the excited- 
state mass is smearing dependent. In Figs. [10] and [TTl 
the ground-state mass reveals no significant dependence 
on smearing. However, the masses of the excited state 
show a distinct systematic dependence on the smearing 
radius. Horizontal lines are drawn at the radius ~ 3.52 
(no. of sweeps — 35) to aid in illustrating the absence of 
source invariance. Fig. [11] presents similar results for the 
lightest quark mass considered in this analysis. 

One might search for an optimal level of smearing 
where the excited state mass plateaus indicating over- 
lap with a neighbouring state is minimized. However, 
there is no evidence of a plateau in Figs. [TUjandfTTI 



This behaviour is also present for our 3x3 correla- 
tion matrix analysis as illustrated in Fig. [T^] Here it is 
found that the two excited states are almost degenerate 
and display a similar dependence on the source smearing 
parameters. 

Thus, we must conclude that the standard analysis of 
the 3x3 correlation matrix of xi, X2 and Xi interpola- 
tors is insufficient to isolate the energy eigenstates. The 
first excited-state mass revealed here is due to a linear 
combination of mass eigenstates and therefore is likely to 
sit high relative to the ffist excited cigcnstate mass. 

In Fig. [TS] we present results for a variational anal- 
ysis of smeared-smeared correlation functions. The re- 
sult from the 3x3 correlation matrix analysis is shown 
in Fig. [13] In this case we also observe two nearly- 
degenerate excited states. While there is some sugges- 
tion of a plateau in this case, it seems unlikely to us that 
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FIG. 5: (Color online). Mass of the nucleon (A''^ ) from the projected correlation function as shown in Eg 1321 fleft) and from 
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to the diagonalization of the correlation matrix for each set of variational parameters tstart (shown in major tick marks) and 
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FIG. 6: (Color online). As in Fig. O but for the source smeared to point sink correlation functions with the number of 
Gaussian smearing sweeps at 7, which corresponds to a smearing radius of 1.5831 in lattice unit. 



the masses revealed here are true eigenstate masses. 



CONCLUSION 



It is important to consider the impact of the finite vol- 
ume of the lattice on our observations. It is well known 
that the eigenstate energies have a volume dependence 
and will change as one changes the volume ,20, 21i|. 
However, as we consider only one fixed volume, the eigen- 
state energies are fixed. Thus the variation of the excited 
state mass revealed as our interpolating fields change can 
only be due to a superposition of eigenstates in the effec- 
tive mass function. 



In this study we have defined and demonstrated a ro- 
bust technique for the analysis of correlation function 
matrices. We have observed that the eigenvectors de- 
scribing the optimal overlap of interpolating fields for 
isolating the first excited state are insensitive to the pa- 
rameters of the eigenvector analysis. This approximate 
invariance of the eigenvectors is in sharp contrast to the 
eigenvalue itself. The latter changes significantly as the 
starting time and the change in time is varied. To create 
a robust technique for the extraction of the excited state 
mass, we exploit the invariance of the eigenvectors and 
construct an eigen-projected correlation function. This 
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FIG. 7: (Color online). Eigenvector values of «f (left) and uf (right), as shown in Eq. (|32p . for the correlation matrix analysis 
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FIG. 8: (Color online). As in Fig. [T] but for the correlation matrix analysis of Fig. |6]and for the excited state. 
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FIG. 9: (Color online). Eigenvector values of v" and u°', for the 3x3 correlation matrix analysis of xi, X2 and X4 interpolators, 
for point source to point sink correlation functions and for the pion mass of 797 MeV. The left figure corresponds to the 
contributions of the interpolators to the ground state, while middle and right figures are for the first and second excited states. 
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FIG. 10: (Color online). Mass of the nucleon (A^2 ) from the projected correlation functions for the pion mass of 797 MeV, 
for point (leftmost point) and for the source smeared to point sink correlation functions (all other points) with rms radii 
0.6897, 1.0459, 1.5831, 2.0639, 2.3792, 3.0284, 3.5237, 4.1868, 5.0067, for 2 x 2 correlation matrices of Xi, X2 (left) and Xi, X4 
(right). Horizontal lines are drawn through the points corresponding to radius 3.5237. 
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FIG. 11: (Color online). As in Fig. [TO] but for the pion mass of 541 MeV (lightest). 
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FIG. 12: (Color online). As in Fig. 1101 but for the 3x3 correlation matrix of xi, X2 and X4 interpolators. The left figure 
corresponds to the pion mass of 797 MeV, whereas the right figure corresponds to a 541 MeV pion mass. 
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lines are drawn through the points corresponding to radius 
3.0284. 



correlation function is analysed using standard analysis 
techniques. 

To reduce human intervention in the fitting procedure, a 
fitting algorithm has been developed which is governed 
by specific fitting criteria based on the maximization of 
the fit window and the minimization of the value of the 
X^/dof while commencing at the earliest time-slice pos- 
sible. 

This study has shown for the first time that the ex- 
cited state masses are fermion-source smearing depen- 
dent for all three types of smearing combinations; i.e., 
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FIG. 14: (Color online). As in Fig. [13 but for the 3 x 3 

correlation matrix of xi, X2 and X4 interpolators. 



the smeared source with point sink, the point source with 
smeared sink and smeared-smeared combinations. This 
is a somewhat unexpected result given that the ground 
state mass is independent of smearing. 
All our 3x3 correlation matrix analyses provide two 
nearly degenerate excited-state masses. However, our 
concern is that these two masses correspond to strongly 
mixed QCD eigenstates that our analysis using the stan- 
dard correlation matrix of xi j X2 and Xi interpolators is 
unable to resolve. 

In particular, this technique has been used by several re- 
search collaborations to determine the mass of the Roper 
resonance [Sll, . Remarkably, mass estimates based on 
these correlation matrix techniques tend to sit high rel- 
ative to other approaches. This investigation provides a 
plausible explanation for these discrepancies. 
Finally, it is clear that changing the smearing level of the 
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fcrmion source and sink changes the relative overlap of 
the superposition of the true eigenstates of QCD. Thus it 
would be interesting to use the robust analysis techniques 
presented here with large correlation matrices built not 
only on the xii Xi ^^id X4 interpolators but also on sev- 
eral levels of fermion-source and -sink smearing. This 
will be the subject of a future investigation. 
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